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Abstract 

A relativistic few-body theory is formulated using point form quan- 
tum mechanics, in which all of the interactions are in the four-momentum 
operator and Lorentz transformations are kinematic. The four-momentum 
operator is written as a product of mass and four-velocity operators, 
where the mass operator is the sum of free and interacting mass oper- 
ators. Interacting mass operators are constructed from vertices, prod- 
ucts of local field operators, evaluated at the space-time point zero. 
Matrix elements of such mass operators, evaluated on four-velocity 
eigenstates of a truncated Fock space, which is the space of the few- 
body theory, are shown to behave like relativistic potentials. Examples 
for a simple vertex are given]^. 

1 Introduction 



One of the goals of few-body nonrelativistic nuclear physics has been to use 
phenomenological potentials determined by two-body systems (for example 
the deuteron and nucleon-nucleon scattering) to predict the behavior of three 
or more body systems. Similarly, current operators determined by one and 
two body systems are used to predict form factors for three or more body 
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systems. In the energy domain where nonrelativistic quantum mechanics is 
vahd this program has been quite successful 

To carry out a similar program at higher energies requires introducing 
some form of relativistic quantum mechanics, not only because of the higher 
energies involved, but also because particle ( for example pi meson) produc- 
tion becomes important. One of the ways to generate a Poincare covariant 
theory is via the so-called Bakamjian-Thomas procedure where inter- 

actions for the few body system are put into a mass operator, which takes 
the place of the Hamiltonian. While this procedure is indeed relativistic, it 
generally violates cluster separability properties, wherein if two clusters of 
interacting particles are widely separated, there should be no interactions 
between the separated clusters. It has been shown that cluster separability 
properties can be satisfied by introducing appropriate unitary operators 
but this procedure is unwieldy, and in any event it is not entirely clear how 
to extend it to deal with particle production. 

Quantum field theory on the other hand provides a setting in which par- 
ticle production occurs naturally, as does cluster separability. However when 
the Fock space on which the quantum field theory is formulated is truncated 
to a finite number of degrees of freedom, relativity is lost. The goal of this 
paper is to construct mass operators acting on a truncated Fock space for an 
n-body system. Up to a given truncation the mass operator should determine 
all the possible reactions that were allowed by the field theory. Further the 
many-body theory should satisfy the following properties: 

-For a given truncation of Fock space the dynamics should arise from a 
many-body mass operator. 

-Generators built out of the mass operator (and other kinematic opera- 
tors) should satisfy the Poincare commutation relations. 

-Forces come from vertex interactions, which generate the mass operator, 
yet should look like relativistic "potentials". 

-The theory should be Lorentz covariant. 

-Bound state problems are solved as eigenvalue-eigenvector problems, 
from which renormalized masses can be calculated. 

The theory to be formulated here is intermediate between a relativistic 
quantum mechanics for a fixed number of particles, and quantum field theory, 
as a theory with infinite numbers of degrees of freedom, and incorporating full 
particle production and vacuum structure. It is reminiscent of earlier models 
such as the Lee model |^ , with its mass renormalization and limited particle 
production. Though the original Lee model was not relativistically covariant. 
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Fuda has shown, using the front form, how to make it properly relativistic. 
The theory developed here differs from this generalized Lee model in that 
it allows for an arbitrary number of produced particles, up to a prespecified 
limit, given by the original truncation. As will be discussed in the conclusion 
there are still open questions about the cluster properties for such a theory. 

The procedure for obtaining a mass operator will be carried out in the 
context of point form relativistic quantum mechanics , in which all of the 
interactions are put in the four-momentum operator and the Lorentz gener- 
ators are kinematic. The four-momentum operator is written as the product 
of a mass operator times a four- velocity operator, which is the Bakamjian- 
Thomas construction in the point form. The interacting part of the mass 
operator is obtained from matrix elements of a vertex, products of local field 
operators evaluated at the space-time point zero, using velocity states, eigen- 
states of the four-velocity operator, in which the matrix elements are diagonal 
in the four-velocity. 

In section 2 the elements of point form quantum mechanics are reviewed 
and both free and interacting mass operators are constructed, using velocity 
states, diagonal in the overall four-velocity of an n particle state. In section 
3 the techniques discussed in section 2 are applied to the simple example of 
a scalar "nucleon", scalar "pion" vertex. 

2 Point Form Relativistic Quantum Mechan- 
ics and Mass Operators 

In point form relativistic quantum mechanics all interactions are put in the 
four-momentum operator. The Lorentz generators contain no interactions, so 
the point form is manifestly Lorentz covariant. In fact the Lorentz generators 
will almost never be used; instead global Lorentz transformations will define 
the transformation properties of operators and states. If ?7a is a unitary op- 
erator representing the Lorentz transformation A, then the four-momentum 
operator must satisfy the following point form equations: 

[P,.PA = (1) 

UaP.U^' = (A-XP^- (2) 

These equations are simply one way of writing the Poincare commutation 
relations in which the relations of the four-momentum operators are empha- 



3 



sized. The mass operator is defined to be M := \JP ■ P and must have a 
nonnegative spectrum. 

Since the four-momentum operators are the generators of space-time 
translations, they can be used to define the relativistic generalization of the 
time dependent Schrodinger equation, 

% - (3) 

where '^x is an element of the Hilbert space. x(= x^) is a space-time point 
and acts as four parameters, playing the same role as the time parameter 
in the nonrelativistic time dependent Schrodinger equation. From Eq.3 it 
follows that \E'a; satisfies a generalized Klein-Gordan equation, 

' ^ ^ +M2)v1/, = 0, (4) 



where M is the mass operator. It is Eqs.3,4 that take the place of equations of 
motion for operators. If has no explicit space-time dependence then Eq.3 
can be written as an eigenvalue equation for the four-momentum operator: 

P^^ = p^^. (5) 

The simplest way in which the point form equations can be satisfied is 
with single particle irreducible representations of the Poincare group. Then 
the Hilbert space is L'^{P?) , where is the dimensional spin 

space for a particle of spin j, and the state |p, a > transforms under Poincare 
transformations as 

P^^\P,(^> = Pt,\p,(^>, (6) 
UA\p,a> = Y.\^p,a' >Dl,^^{Rw), (7) 

Rw : = B-\Av)AB{v). (8) 

Here the four-momentum p satisfies the mass constraint p ■ p = m? and 
the four-velocity v is defined by w = p/m. B{y) is a boost (see reference 
0), a Lorentz transformation that takes the rest four- momentum, p^'^*** = 
(m, 0, 0, 0) to the four-momentum p: B{y)p''''^'^* = p. Rw is a Wigner rotation 
and -D-'O is a matrix element of the rotation group. It is clear that the 
operators Eqs. 6,7 satisfy, by construction, the point form equations. 
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Many-particle states with the same transformation properties as the single 
particle ones are conveniently obtained by introducing creation and annihi- 
lation operators. Let a^(p, cr) be the operator that creates the state \p,a > 
from the vacuum. If a{p, a) is its adjoint, these operators must satisfy the 
following relations: 

[a(p,a),at(p',a')]± = E5^p-p)6^y (9) 
C/„at(p,a)t/-i = e*PV(p,a), (10) 

P.ifr) = T.J^P.^Hp,<^Hp,<t), (11) 

UAa^{p,a)Uj[' = Y.'^\Ap,a)Dl,^^{Rw). (12) 

Here P^{fr) is the free four-momentum operator and plays a role analogous 
to the free Hamiltonian in nonrelativistic quantum mechanics. Again it is 
straightforward to show that satisfies the point form equations. Ua in 
Eq.lO is the unitary operator representing the four-translation a. 

To prepare for the construction of interacting mass operators, it is con- 
venient to introduce velocity states, states with simple Lorentz transforma- 
tion properties. If a Lorentz transformation is applied to a many-particle 
state, (Ti...p„, cr„ >= a\pi,ai)...a\pm(Jn)\Q >, then it is not possible to 
couple all the momenta and spins together to form spin or orbital angular 
momentum states, because the Wigner rotations for each momentum state 
are different. However, velocity states, defined as n-particle states in their 
overall rest frame boosted to a four-velocity v will have the desired Lorentz 
transformation properties: 

\v,ki,IJ,i>: = UB{v)\ki,IIi...kn,IIn> (13) 

= ^|pi,ai...p„,a„>n^i;,^,(^w^J- (14) 

— * 

UA\v,ki,IIi> = UAUB{v)\kl,IIl...kn,IIn> 

= UB{Av)URy^,\ki, IJ,i...kn, IJ'n > 

= J2\^v,Rwh,l^'i>IlD'}ARw). (15) 

Unlike the Lorentz transformation of an n-particle state, where all the Wigner 
rotations of the D functions are different, in Eq.l5 it is seen that the Wigner 
rotations in the D functions are all the same and given by Eq.8. Moreover the 
same Wigner rotation also multiplies the internal momentum vectors, which 
means that for velocity states, spin and orbital angular momentum can be 
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coupled together exactly as is done nonrelativistically (for more details see 
reference 0). The relationship between single particle and internal momenta 
is given by pi = B{v)ki,J2 h = 0; Eq.l4 provides the link between velocity 
and many-particle states, where is given by replacing p by ki and A by 
B{v) in Eq.8. From the definition of velocity states it then follows that 

Vf,\v,ki,fii> = Vf,\v,ki,fii>, (16) 
M{fr)\v,ki,fii> = mn\v,ki,fii>, (17) 
P^{fr)\v,ki,Hi> = mnv''\v,ki,jji>, (18) 



with m„ = J2 \l f^f + ki the 'mass' of the n-particle velocity state and 
Pfiifr) = M{fr)Vf^. On velocity states the free four-momentum operator 
has been written as the product of the four-velocity operator times the free 
mass operator, which is the Bakamjian- Thomas construction for the point 
form. More generally four- momentum operators are written as = MV^, 
where the four- velocity operator is defined by := ^-j^j-^', is diagonal on 
velocity states, as seen in Eq.l6. The mass operator is the sum of free and 
interacting mass operators, M = Mf^ + Mj; if the mass operator commutes 
with the four-velocity operator and Lorentz transformations, then the point 
form Eqs.1,2 will be satisfied. 

The goal now is to construct an interacting mass operator on a truncated 
Fock space, where the Fock space is formed from creation and annihilation 
operators. The truncated Fock space is the space of the many-body system, 
a system with a finite number of degrees of freedom. For example if the 
creation and annihilation operators are those for (bare) pions and nucleons, 
the truncated Fock space could be chosen as the direct sum of two nucleon 
spaces with zero to n pions, appropriate for the scattering of two nucleons as 
well as the production of pi mesons. 

The same creation and annihilation operators can also be used to form 
local fields; the general procedure for this construction is given in reference 
H and is reviewed in detail in one of the papers in this series, dealing with 



constructing current operators for arbitrary spin particles [^. Moreover 
reference shows how to take products of such arbitrary spin fields and 
couple them together to form Lorentz scalar densities, which, following this 
reference will be denoted by T-C{x). Examples of such scalar densities of 
products of local fields are the familiar ^^^"$(1) for pion-nucleon vertices and 
^'j^'^A^ for electron-photon vertices. 
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The fact that 7i{x) is a Lorentz scalar density means that U\H{x)U]^^ — 
T-C{Ax). Therefore at the space-time point a; = 0, UaH{0)UJ^^ = 7Y(0), so 
that 7i(0) is a Lorentz scalar. The velocity state matrix element of 7i(0) 
on the truncated Fock space, < v' kliJ,[\T-C{0)\vkiiJ,i > where the initial 

and final four-velocities are the same, can be viewed as the kernel of an 
interacting mass operator on the truncated Fock space, because it commutes 
with the four-velocity operator (since it is diagonal in the four-velocity, by 
construction) and commutes with Lorentz transformations. 

However, since 'H(O) is constructed from products of fields, its velocity 
state matrix elements may not give kernels of well defined operators on the 
truncated Fock space. In order that the interacting mass operator be well 
defined on the truncated Fock space, matrix elements of ?i(0) are multiplied 
by "form factors", functions of {p — pY = {m v — mvY, where m and m 
are the free masses of the many-particle state, given in Eq.l7. But since 
the initial and final four-velocities are the same, the form factor will be a 
function of the magnitude of mass differences only, namely Am := |m — m|. 
The definition of the interacting many-body mass operator is then given by 

Ml-. = <v',klfi[\n{0)\v,k„fi,>l>^J{Am), (19) 

and by virtue of the form factor /(Am) is a well defined operator on the 
truncated Fock space. It commutes with Lorentz transformations and the 
four- velocity operator, so that = {Mj^ + Mi)V^ satisfies the point form 
equations, Eq.1,2 and provides the starting point for a relativistic many-body 
theory. 

It should be noted that Mj actually is independent of the four- velocity: 

Mi = < v',fc;,/i;|7f(0)|t;,fc„/i, > |„.^J(Am) 

= < k!-^^J!^...\UJ^l^^n{())UB(v)\kl^ll... > |„'=j(Am) 

= < k[i,[...\nmk,f,,... > |„'=J(Am), (20) 

where the fcj and k[ are the internal momenta defined in Eq.l3 and use has 
been made of the fact that 7i(0) is a Lorentz scalar. Also, since under 
Lorentz transformations the internal momenta are rotated by the Wigner 
rotation of Eq.l5 the mass operator is rotationally invariant and satisfies 
all the properties of a "potential". Finally, since Mj arises from products 
of fields, it always has only off-diagonal matrix elements in the truncated 
Hilbert space. 
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3 An Example: The Scalar Nucleon-Pion Ver- 
tex 



To show how mass operators and generahzed Lippmann-Schwinger equa- 
tions result from vertices, this section works out the example of a scalar 
"nucleon" interacting with a scalar "pion" . The vertex is of the form H{x) = 
g'i>\x)'^{x)(l){x) where g is a coupling constant. The interacting mass oper- 
ator can be written from Eq.l9 as 

Mi = gf{Am)<v,h,^^\'i^\0)^{0)<j>{Q)\v,k■,|J,[>, 

and various baryon number sectors evaluated. 

As a first truncation consider the baryon number one sector, with the 
Hilbert space the direct sum of bare nucleon plus bare nucleon and bare 
pion. Then the total mass operator can be written as a two by two matrix 
operator in the direct sum space: 



M = 



mAr(O) gK^ 
gK D^iO) 



(21) 



where the internal momentum variables are defined by 

Pn — B(v){ujN{0),k) and p-,^ — B{v){uJt,, —k) with a;Ar(0) = J m^(0) -|- k? 



and cUtt = v rnl ^k"^. rriN^Q) is the bare nucleon mass and 



1^2(0) := yfn%{0) + k'^ + yml. + k'^ = tUAr(O) -|- a;^, a diagonal operator in 
the two particle space, is the relativistic two particle mass. The interacting 
mass operator that connects the two particle space to the one nucleon space 
is the kernel 

< k\K\ > = /(Am) 

= f{ujMiO)+uj^-m^{0)). (22) 

On the direct sum space the problem to be solved is the mass operator 
eigenvalue-eigenvector problem, Mcf) — m0, where, for the bound state prob- 
lem, m = is the physical nucleon mass, while for the scattering problem, 
m = is the invariant relativistic energy for pion-nucleon scattering. Des- 
ignating the components of the physical one nucleon state as (01^,02^), the 
mass operator eigenvalue equation becomes 

mAr(OX + ^(i^,0^) = mAT^f (23) 
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N 



N 



0: 

g\K,(mN-D2{0))-'K) = 



^(m^ - L'2(O))-'X0f 
mAT - mAr(O), 



(24) 
(25) 
(26) 



where in the last equation, 0f has been cancelled from both sides of Eq.26 
because it is a number (to be determined by normalization requirements). 
For a given truncation, Eq.26 provides the relationship between the bare and 
physical nucleon mass. If m^r is taken to be the physical nucleon mass, then, 
for a given form factor, the bare nucleon mass, mjv(O), is fixed. Written out 
more explictly, Eq.26 becomes 



d^k |/(a;iv(0)+a;^-miv(0))p 



= mN-mN{0), (27) 



and / is used to regularize the large momentum components in the integral. 
The nucleon wave function is normalized to one; that is, |0]^p+(0^, 0^) = 1, 



where 0^ is given in Eq.25. 



To formulate the pion-nucleon scattering problem, the mass operator, 
Eq.21, is written as the sum of the free mass operator, Mq, and the interacting 
mass operator, M/, M — Mq + Mj, where 



Mo = 
Mj = 



miv(O) 

^2(0) 

gK^ " 
gK 



(28) 
(29) 



The relativistic Lippmann-Schwingcr equation can then be written as an 
integral equation on the direct sum space with scattering wave function ^' = 
(*i,^'2): 



^ = $ + (Vs-Mo + ie)"^M7* 
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^/s - mAr(O) 



02 + 



- mAr(O) 



{^/^-D2{Q))-'K{K, *2), 



(30) 
(31) 

(32) 



where the second term in Eq.32 looks like a separable potential with a cou- 
pling constant g^j^^fs — m7v(0)). That is, in this simple truncation the 
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Lippmann-Schwinger equation can be solved exactly to compute the pion- 
nucleon scattering amplitude. $ = (0,02)) where 02 is the free pion-nucleon 
momentum state, and all ^/s factors in Eqs. 30,31 should be replaced by 
^/s + ie. 

The function f{Am) used to relate bare and physical mass, Eq.27 and 
generate the separable potential, Eq.32, can also be used in other baryon 

number sectors. For baryon number two, consider the truncated space con- 
sisting of two bare nucleons plus two bare nucleons and bare pion, a direct 
sum of two particle and three particle spaces. In this case the mass operator 
has the form 



M 



D2{0) gK^ 
gK Ds{0) 



(33) 



where the internal momenta are related to single particle momenta by the fol- 
lowing relations: p'l = B{v){cu, k),p'2 = B{v){cu, —k), with a; = ^m^(O) + fc^; 



Pi = B{v){ui,ki),p2 = B{v){uj2,k2),p^ = B{v){uj^,k^), with uji = ^Jm%{0) + kf 
and ki + k2 + k^^ = 0. The diagonal operator D2{0) = 2uj is the bare two 
particle relativistic mass, while the diagonal operator 1^3(0) = Ui + uj2 + uJtt 
is the bare three particle mass. The two particle to three particle transition 
kernel is given by 

< kik2\K\k > = /(c<Ji + CJ2 + i^TT — 2c<j) 

{5%k + h) + S^k - h) + 5{k + 4) + 5{k - 4)),(34) 

where the argument in / is (the magnitude of) the mass difference between 
the three and two particle states. The four delta functions are a consequence 
of the two nucleons being identical particles. 

With the mass operator given in Eq.33 a "deuteron" bound state can be 
calculated as a two component mass operator eigenvector, corresponding to 
the deuteron mass niD'- 



M0^ 




(35) 


D2m2+gK^<t>^ 


= mD(t)2 


(36) 


gK<P^ + Dsmf 


= rnD(j)3 


(37) 




= giniD - 


Ds{Q))-'K(j)^ (38) 


D2(O)0f + g^K^imo - D^y^KcP^ 


= mD(p2- 


(39) 
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Though the first term in Eq.39 is the bare relativistic kinetic energy, the 
second term, which looks hke a potential energy term, contains the unknown 
deuteron mass eigenvalue; because the component of the deuteron wavefunc- 
tion in the three particle sector has been eliminated, Eq.39 is not of the usual 
form of an eigenvalue equation, but is known in coupled channel problems. 
The deuteron wave function is normalized by writing ||0f |p + ||0f |p = 1. 

The mass operator, Eq.33 can be used to calculate scattering in the 
nucleon-nucleon and nucleon-nucleon-pion channels by again writing the mass 
operator as a sum of the free mass operator plus interacting mass operator, 
M = Mo + Mj: 



Mo = 

Mi - 

* = 

^2 = 

*3 = 

*2 = 



1^2(0) 
1)3(0) 

gK^ " 
gK 

$ + ( - Mo + ie) 



(.2 + - D2m-'gK^-^s: 
-Dsm-'gK^2 
+ g\^s - D,{0))-'K^{^s D-Mr'K^ 



2, 



(40) 

(41) 

(42) 
(43) 
(44) 
(45) 



where $ = (02,0), and the driving term is the two free nucleon momentum 
state, 02 = S^{k — kiniuai)- Again Eq.45 has the form of a driving term plus 
"potential", from which it is possible to calculate scattering amplitudes for 
the two processes mentioned above. 

Finally, the baryon zero sector gives the ground state and physical pion 
mass. However in the truncation where there are no nucleon-antinucleon 
pairs, and only the direct sum of Fock vacuum, and one or more bare pion 
states, there is no pion mass renormalization and the physical vacuum is the 
Fock vacuum. If however a nucleon-antinucleon two particle space is added 
to the Fock vacuum along with a nucleon-antinucleon-pion space and a one 
and two pion space, the mass operator will have off-diagonal terms linking 
the Fock vacuum to the nucleon-antinucleon-pion Hilbert space, and pion 
mass renormalization will occur, as well as pion-pion scattering. 

More generally for a given truncation the mass operator will be a quint- 
diagonal matrix operator for which bound and scattering states can be cal- 
culated. If the function / occurring in the off-diagonal kernal is fixed, then 
for the given truncation all channels will be determined by /. 
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4 Conclusion 



This paper has shown how to compute interacting mass operators from a 
given vertex, by muhiplying velocity state matrix elements of the vertex op- 
erator by a form factor that makes the mass operator a well defined operator 
on the truncated Fock space appropriate for a few-body system. The total 
mass operator, the sum of free and interacting mass operators, has a dis- 
crete spectrum that gives the bound states, and continuous spectrum the 
scattering states. And just as the total momentum and spin projection are 
extracted from an isolated many-body nonrelativistic system, so in the point 
form the overall four-velocity and spin projection are extracted from an iso- 
lated many-body relativistic system. Since the interacting mass operator 
comes from a vertex operator, the interactions are all due to the exchange 
of particles. Yet, as shown in section 3 the coupled channel problem can 
be rewritten so that the relativistic Lippmann-Schwinger equations contain 
relativistic "potentials". 

There are a number of applications of the formalism discussed here. The 
most obvious is a more realistic treatment of the pion-nucleon system, in 
which the spins are correctly taken into account. Then the interacting mass 
operator will contain spinors, but, as can be seen from the way in which fields 
for arbitrary spin particles can be coupled together to form Lorentz scalar 
densities [^, the coupling can also include the spin 3/2 delta resonance. Lee 
and coworkers |]lT|, as well as Fuda |jl2[ have written T matrix equations 



for pion-nucleon scattering, using basically the instant form of relativistic 
quantum mechanics, in which vertices are used to generate interactions in 
the Hamiltonian. A similar procedure should be possible in the point form. 

A second application is to the study of the meson and baryon spectrum, 
arising as bound states of quarks. While the work of Plessas et al [T^ can 
be interpreted as a point form mass operator, the hyperfine part of the mass 
operator only indirectly comes from a vertex. Closer in spirit to what is devel- 
oped here is the dissertation of Krassnigg |]T^ who studies the vector mesons 
as bound states of quark-antiquark pairs. The truncated Hilbert space is 
the direct sum of quark-antiquark and quark-antiquark, pseudoscalar me- 
son spaces, and the hyperfine interaction is given from a vertex of the form 
1-L{x) = 'l'(x)75\l/(x)0(x), where ^(x) is the quark field and 0(x) the pseu- 
doscalar meson field. The main difference with the mass operator discussed 
in this work is that the free mass operator is replaced by a harmonic os- 
cillator mass operator to simulate quark confinement. But otherwise the 
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coupled channel equations are very similar to those of section 3, and result, 
for example, in excited states having widths and shifted masses. 

For a given truncation there is a definite mass renormalization, which 
however will change if the truncation changes. Already at this point viola- 
tions of cluster properties will occur, for if a set of bare particles are moved 
far away from another set, the mass renormalization for the separated set will 
not be the same as with all the particles. But for a given truncation and mass 
renormalization, it should be possible to construct unitary operators (pack- 
ing operators,see reference 0) that give the desired cluster properties for the 
many-body system. Further there is reason to believe that the packing oper- 
ators for vertex interactions should be relatively straightforward to calculate. 
The desirable cluster properties come about because / has as an argument 
the difference of the relativistic masses of the initial and final sets of parti- 
cles. The mass operator is the product of this function multiplied by velocity 
state matrix elements of the vertex operator. For a given vertex only three 
particles interact while the momentum of the others are unchanged, resulting 
in the cancellation of energies in the argument of /. Then the argument of 
/ depends only on the momentum of the interacting particles and not the 
others. Because of this peculiar feature of the argument of / in the point 
form, the cluster properties may be relatively simple. But this remains to be 
investigated. 

More generally as pointed out in the introduction, the coupled channel 
mass operators that result from vertex interactions allow for a limited amount 
of particle production, depending on the truncation chosen. Such theories 
are intermediate between relativistic direct interaction theories for a fixed 
number of particles and quantum field theories with infinite degrees of free- 
dom. Moreover, quantum field theories satisfy cluster separability properties 
(see for example, 0), while relativistic direct interaction theories for a fixed 
number of particles require packing operators to satisfy cluster properties 
1^. While the relativistic coupled channel theory discussed in this paper 
has limited particle production and requires mass renormalization, the way 
mass renormalization depends on cluster properties remains an open ques- 



tion. This issue is particularly relevant in light of a theorem by Aks |15 
that says in a four dimensional quantum field theory particle production and 
vacuum polarization necessarily go together. 

The procedures given in this paper for constructing strongly interacting 
mass operators from vertices can also be applied to electromagnetic vertices, 
where ?i(0) is of the form J^{0)A^{0)] in the following papers current oper- 
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ators are constructed for particles with arbitrary spin and form factors |TD 



as well as analyzing the Gupta-Bleuler formalism appropriate for the above 
coupling |]T6[. Finally the last paper in this series [|I^ shows how the free cur- 
rent operator is modified from its one-body form to generate the appropriate 
many-body current operator in the presence of interactions. The goal then, 
in this series of papers is to see whether the nonrelativistic methods that have 
been so successful in low energy nuclear physics, where the Hamiltonian is 
the sum of matter, photon, and electromagnetic coupling Hamiltonians, can 
be generalized to higher energies where relativity is required and the mass 
operator is similarly the sum of matter, photon and electromagnetic coupling 
mass operators. 



5 Appendix 

In this appendix the interacting mass operator is obtained from 7Y(x) by inte- 
grating over the forward hyperboloid to get the interacting four-momentum 
operator. Let 7i(x) be a polynomial in free fields which is a Lorentz scalar 
density, UaH^x)!!^^ = H^Ax). Then an interacting four-momentum oper- 
ator which satisfies the point form equations can be defined by integrating 
T-C{x) over the forward hyperboloid: 

P^(/) : = J d^x6{x ■ X - T^)Xf,e{xo)n{xy, (46) 
[P^(/),P.(/)] = 0. (47) 

The components of the interacting four-momentum operator commute with 
one another since the commutator of ?i(x) with 7^(?/) is zero if the space- 
time points X and y are space-like separated. But both these points lie on the 
same forward hyperboloid specified by r and hence are space-like separated. 
It also follows from Eq.46 that P^{I) transforms as a Lorentz four- vector, 
since 7i(a;) transforms as a Lorentz scalar density. Hence the interacting 
four-momentum operator satisfies the point form equations, Eq.1,2. 

The total hadronic four-momentum operator is the sum of free and in- 
teracting four-momentum operators, P^(/i) = P^(/r) + Pfi{I)- Again it can 
be shown that the total four-momentum operator satisfies the point form 
equations. The Lorentz transformation part follows from the fact that both 
terms transform as four-vectors. To show that the components of the total 
four-momentum operator commute, it suffices to note that the components 
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of the free part commute among themselves (from Eq.ll), as do the compo- 
nents of the interacting part (Eq.47). Thus what must be shown is that the 
sum of the cross term commutators give zero. Now from the definition of 
7i{x) as a polynomial in free fields, it follows that Ua'H{x)U~^ = 7i{x + a). 
But Ua is the exponential of P{fr) ■ a, so if a is made infinitesimal, then 

[P,(/r),P^(7)] = ld^xS{x-x-ry{xo)x^-^n{xy, (48) 
+[P^(/),P.(/r)] = [P,(/r),P.(/)]-[P.(/r),P^(7)] 

/d 
d!^x5{x ■ X - t'^)9{xo){x^— - x^—)H{x) 

= 0, (49) 

since 7i(x) is a Lorentz scalar density. Thus the total four-momentum oper- 
ator satisfies the point form equations. 

Given the total four-momentum operator, the sum of Eqs. 11,46, bound 
and scattering states could in principle be found by solving the eigenvalue 
equation, Eq.5. But the interacting four-momentum operator defined in 
Eq.46 is not a well defined operator because it involves the product of local 
field operators at the same space-time point. Moreover, if a restriction to 
a finite degree of freedom system is made, the theory is no longer properly 
relativistic because Eq.47 is not valid on a truncated space. Finally, the 
interacting four-momentum operator given in Eq.47 is never of Bakamjian- 
Thomas form; that is, it can never be written as the product of a mass 
operator times a four-velocity operator. Nevertheless, a Bakamjian- Thomas 
interacting mass operator can be constructed from Eq.46 by examining ma- 
trix elements of the interacting four-momentum operator, where the states 
are velocity states with the same initial and final four-velocities: 

<v,ki,Hi\P^,{I)\v,k\,^^^> = < f;, //i|'H(0)|^;, > 

j d'^xSix ■ X - r'^)x^e{xo)e'^^"'>-'', (50) 

where Am is m — m' and m (respectively m) is the mass of the velocity 
state, as given in Eq.l7. {Am)v is always a timelike four- vector, although 
the sign of the time component may be positive or negative, depending on 
the sign of Am. 
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The integral in Eq.50 is a special function and is evaluated in the following 
paragraph. The relevant point is that, if the initial and final four- velocities 
are arbitrary, the integral is not diagonal in the four- velocity, which is another 
way of saying that an interacting four-momentum operator in quantum field 
theory is never of Bakamjian-Thomas form. 

But Eq.50 can be used to construct an interacting mass operator, gener- 
ated from vertex interactions, by defining 

Mi = <vkifii\n{0)\vk[fi[> f{Am), (51) 

where / is the function defined in Eq.50 and evaluated in the following para- 
graph. If it is replaced with an arbitrary form factor, the mass operator is the 
same as that given in Eq.50, which was the starting point for a many-body 
theory. 

The integral in Eq.50 can be computed as follows. Define the Lorentz 
invariant function as 

J±(p) = j d^x5{x ■ X - T'')e{xo)e^"'"' (52) 
I${Kp) = I^p), (53) 
I^sp) = s-'ltAp) (54) 

From reference ||18[ the function I^{p) is given by 

J+(p) = 7T^r{p-py'/^[N,{T,/^)-ie{po)Mry^)] (55) 
Irip) = ilHpW, (56) 

for p timelike. 

Then for Am > 0, 

j d^x5{x-x-T^)x^e{xQ)e'^^'^^'''' = {iAm)-^d/dv^'I+{Amv) 

= v'^friAm), (57) 

where /,-( Am) = —-K'^T{Am)^'^a~^d/da[a'^Ji{AmTa)+ia~^Ni{AmTa)]a=i 
with a := y/v ■ v. For Am < 0, the desired integral is 

J d^x6{x-x-T^)x^e{xo)e-'^^"'^''-'' = ilAml-^d/dv^I^^^^^^^iv) 

= v'^miAml). (58) 
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The integrals in Eqs. 57,58 are complex conjugates of one another, with both 
depending on the magnitude of the difference between final and initial masses. 
Ji and A^^i are various forms of Bessel functions, as given in reference . 
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